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Анотація. При обробці та передачі інформації у вигляді зображень, важливою є проблема 
зменшення їх спотворення за рахунок різноманітних шумів. Шум знижує якість зображення 
і тим самим сприйняття інформації, що міститься у ньому. Через ці проблеми погіршується 
можливість оцінки інформації, яку можна отримати в результаті аналізу за допомогою як 
візуального так 1 комп'ютерного методів. Процедурою зменшення шумів на зображеннях 
займається область обробки зображень, яка отримала назву відновлення. Не дивлячись на 
перетин цієї області з покращенням зображень, слід відмітити, що остання є у більшій степені 
суб'єктивною процедурою ,у той час як процес відновлення має в основному об'єктивний 
характер. При відновленні робиться спроба реконструювати або відтворити спотворене 
зображення, використовуючи для цього апріорну інформацію про явище, яке викликало його 
погіршення. Методи відновлення грунтуються на моделюванні процесів спотворення і 
використанні зворотних процедур для відновлення початкового зображення. У даній роботі 
запропоновано новий метод оцінки відновлення спотворених зображень, здійснених за 
допомогою просторових i частотних фільтрів. Він полягає у переведенні спотворених 
зображень перед фільтрацією у градації сірого 1 отриманні їх гістограм. Отримана гістограма 
екстраполюється гаусівською кривою і визначається величина середньоквадратичного 
відхилення. Аналогічна процедура здійснюється 1 для відновленого зображення. Аналіз 
проводився на мові програмування Python із використанням бібліотек для роботи 3 
зображеннями Pillow та OpenCV. . Запропоновано параметр оцінки відновлення К , який являє 
собою відношення різниці середньоквадратичних відхилень спотвореного і відновленого 
зображень до середньоквадратичного відхилення спотвореного зображення. Оцінка параметру 
К була здійснена для середньогеометричного фільтру, фільтру серединної точки та фільтру 
Вінера 1 виявилась рівною 0,2; 0,3 1 0,43 відповідно, що корелює з візуальними 
спостереженнями. 

Ключові слова. відновлення зображень,  гаусівський шум, фільтр Вінера, 
середньогеометричний фільтр, фільтр серединної точки, Python. 
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1. Beryn 


При обробці та передачі інформації у вигляді зображень, важливою € проблема 
зменшення їх спотворення за рахунок різноманітних шумів. Шум знижує якість зображення 
і тим самим сприйняття інформації, що міститься у ньому. Через ці проблеми погіршується 
можливість оцінки інформації, яку можна отримати в результаті аналізу за допомогою як 
візуального так і комп'ютерного методів. Процедурою зменшення шумів на зображеннях 
займається область обробки зображень, яка отримала назву відновлення. Не дивлячись на 
перетин цієї області з покращенням зображень 1-8), слід відмітити, що остання є у більшій 
степені суб'єктивною процедурою ,у той час як процес відновлення має в основному 
об'єктивний характер. При відновленні робиться спроба реконструювати або відтворити 
спотворене зображення, використовуючи для цього апріорну інформацію про явище, яке 
викликало його погіршення. Методи відновлення грунтуються на моделюванні процесів 
спотворення і використанні зворотних процедур для відновлення початкового зображення. 


2. Об'єкт і предмет дослідження 


В останні роки це процес обробки спотворених зображень здійснюється за допомогою 
високорівневих мов Python ra СЯ. Мова Python використовується у зв'язку з тим, що вона 
володіє широким спектром пакетів, які надаються для обробки зображень. Її бібліотеки та 
інструменти дозволяють дуже ефективно виконувати завдання у цій галузі. Це зокрема 
стосується бібліотек для роботи з зображеннями Pillow та OpenCV. Перевагою мови СЯ є 
доволі великий функціонал , за допомогою якого можна реалізувати програми з ефективними 
способами обробки зображень. Однак при цьому вона має i недолік - відсутність вбудованих 
методів обробки зображень. Тому програму приходиться створювати практично з нуля, що у 
свою чергу зумовлює нагромадження у коді. У зв'язку з цим, у даній роботі нанесення та 
зняття шуму здійснювалось на мові Python. 


3. Мета та задачі дослідження 


Не зважаючи на велику кількість теоретичних і експериментальних робіт в області 
відновлення зображень, існує ще ряд питань, які вимагають свого вивчення та реалізації. Це, 
зокрема, стосується оцінки ефективності використання тих чи інших методів відновлення, які 
давали би оптимальний результат при мінімальних обчислювальних потужностях а також 
створення надійних комп'ютерних програм для відновлення зображень. Тому, у даній роботі 
розглянуто спотворення зображень гаусівським шумом та зняття цих шумів за допомогою 
деяких основних фільтрів (середньогеометричного, серединної точки та Вінера) а також 
пошук критеріїв оцінки ефективності цих методів. 


4. Аналіз літератури 


Поява шуму на зображеннях зумовлена оцифруванням сигналу, тобто перетворенням 
аналогового сигналу в цифровий та його передачею[ 1-3). Під час процесу передачі сигналу 
інформація може спотворюватися через завади що завжди притаманні каналу зв'язку. Також 
можливі спотворення через атмосферні впливи під час бездротової передачі інформації. До 
основних параметрів шуму відносяться[4-8]: просторова характеристика, наявність кореляції 
між шумом та сигналом та частотні характеристики (Фур'є спектр шуму). На рис. 
представлено модель процесу спотворення/відновлення зображення [1]. Згідно цієї моделі, на 
початкове зображення f(x,y) діє певне спотворення , яке описується оператором Н. Після 
накладання на зображення адитивного шуму, отримується спотворене зображення g(x,y). 


^ 


Завдання відновлення полягає у побудові певного наближення  f(x,y), яке відтворює 
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початкове зображення на основі інформації про спотворене зображення g (x, y), яку описують 
оператором Н, та відомостей про адитивний шум N(x, y). 


ED 


Відновлюючий 
оператор 


(фільтр) 


Спотворюючий 
оператор 


ГО» у) 


СПОТВОРЕННЯ 


ВІДНОВЛЕННЯ 
Рисунок 1. Модель процесу спотворення/відновлення зображення 


Необхідно щоби отримане таким чином наближення якнайкраще співпадало з 
оригінальним зображенням. Спотворене зображення можна представити виразом (1) [1-3]: 


9(х, у) = h(x, у) * Р(х,у) + п(х, у), (1) 


де h(x, y) - функція, що представляє собою оператор шуму у просторовій області 
символ ж - використовується для позначення згортки. 


Оскільки згортка у просторовій області є еквівалентною множенню у частотній області, 
то для частотної області напишемо вираз (2): 


С(и, v) = Н(и, 0)Е(и, v) + N(uv). (2) 


V подальшому будемо вважати, що H є тотожним оператором i маємо справу з 
спотвореннями, які викликані лише наявністю адитивного шуму. До найбільш поширених 
видів шумів відносяться: гаусівський; дробовий (шум Пуассона), шум типу «сіль і перець» 
(імпульсний шум) та спекл-шум [1-4]. 


5. Методи досліджень 


У роботі було використане нанесення гаусівського шуму| 1-3]. Математична простота що 
характерна для роботи з моделями гаусівського шуму (який ще називають нормальним 
шумом) як у просторовій так і частотній області зумовила широке використання цих моделей 
на практиці. Функція густини розподілу імовірностей гаусівської величини 7 задається 
виразом (3): 


_ а-и)? 


202 : (3) 


1 
p(z) = Jomm C 


де Z представляє собою значення яскравості; 

и - середнє значеня випадкової величини Z (математичне очікування), 6- ii 
середньоквадратичне відхилення. Квадрат середньоквадратичного відхилення ще називають 
дисперсією величини 7. Графік цієї функції показаний на рисунку 2. 
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Розподілення Гауса 


u- и utc Z 
Рисунок 2. Графік функції розподілу густини ймовірності функції Гауса 


Якщо густина розподілу випадкової величини Z описується функцією (1), то приблизно 
70% попадають у діапазон (и - 6), (и + 6), i приблизно 90%- у діапазон (и — 26), (и + 26). 


У даній роботі для обробки було використане оригінальне зображення [9] , яке 
представлено на рис.1а у кольорі та градаціях сірого 16 відповідно. Тут також представлена 
гістограма чорно-білого зображення (у відносних одиницях |в.о. |). 


B 


Рисунок 3. Використане для опрацювання оригінальне кольорове зображення (а), 
приведення його до градацій сірого (6) 1 гістограма чорно-білого зображення(в) 


У дальшому усі зображення представляються у форматі jpg. Попередні дослідження 
показали, що зміна формату початкових зображень практично не впливає на отримані 
результати, тому у при обробці зображень використовувався формат jpg. Як видно з puc.36 
гістограма початкового зображення має чітко виражений гаусівський характер. Тому, до неї 
можна використати формулу (1) і визначити середньоквадратичне відхилення (СКВ), яке 
виявилось рівним 6 - 0,4. І хоча використання тут гаусівської моделі оправдано лише 
частково, тим не менше спостерігається якісне співпадіння. 
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Рисунок 4. Спотворене гаусівським шумом початкове чорно-біле зображення (а) та його 
гістограма (6) 


На рис4 показано спотворене гаусівським шумом початкове чорно-біле зображення (а) ra 
його гістограма(б). Як видно з порівняння рисунків 3 1 4, гаусівський шум суттєво спотворив 
початкове зображення. Це відобразилось 1 на гістограмах. При цьому величина 
середньоквадратичного відхилення у порівнянні з початковим зображенням збільшилась від 
0,4 до 1,0 (у в.о.) а математичне очікування змістилось вправо до величини 0,25. 

Відновлення за допомогою середньогеометричного фільтру здійснювалось 1з 
використанням виразу (4) [10]: 


f(x, y) = [Поов асо |" , (4) 


де Syy означає прямокутний окіл (множину координат точок зображення ) розміром MXN 
із центром у точці (x, y). 


Процедура фільтрації передбачає обчислення  середньогеометричного значення 
спотвореного зображення g(x,y) в околі Syy. Значення відновленого зображення у кожній 
точці (х,у) є коренем степені MN з добутків значень у точках околу Syy. Як правило, для 
посилення ефекту фільтрації використовують вираз (5): 


f(x,y) = exp [Е 2 stes, In g(s, Э] (5) 


Тут до значень яскравостей спотвореного зображення застосовується логарифмування та 
сумування з усередненням , після чого застосовується функція експоненти. 

Використання фільтру серединної точки полягає у обчисленні середнього значення між 
максимальним 1 мінімальним значенням у відповідному околі, як наведено виразом (6) [11]: 


^ m i я 
70) = 5. max „9 С НА (6) 


Такий фільтр об'єднує у собі методи порядкових статистик 1 усереднення. Він найкраще 
працює при наявності таких випадково розподілених шумів як гаусівський або рівномірний. 

Фільтр Вінера поєднує у собі врахування властивостей спотворюючої функції 1 
статистичних властивостей шуму у процесі відновлення зображень[13, 14,15]. Метод 
грунтується на розгляді зображень 1 шумів як випадкових процесів і завдання ставиться 
наступним чином: знайти таку оцінку Ё для неспотвореного зображення /, щоби 
середньоквадратичне відхилення цих величин одна від одної були мінімальними. 
Середньоквадратичне відхилення задається виразом (7): 
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e -go-5. D 


де ЕС) означає математичне очікування аргументу. 


Допускається, що виконуються наступні умови: шум і неспотворене зображення не 
корелюють між собою; або шум або неспотворене зображення мають нульове середнє 
значення; оцінка лінійно залежить від спотвореного зображення. При виконанні цих умов, 
мінімум середньоквадратичного відхилення (5) досягається на функції, яка задається у 
частотній області за допомогою виразу (8) [13]: 


Н"(и, 2)5; (и, v) 


Кит є S Га x)|H Qu, v)? + S, Cu, v) 


G(u,v) = 


H*(u,v) 
S (и, v) 
5 (и, v) 


G(u,v) = 
|H (u, v)|? + 


1 |H* Quo) 


H(uv) 251007) 
ІН(и,»)| MIO 


G(u, v) , (8) 


де H(u, v) - спотворюючи функція (її частотне представлення ); 

Н" (и, v) - комплексне спряження спотворюючої функції; 

|H (u, v)|? = H*(u, 2)Н(и, v); 5, (и, v) = |N (u, v)|? - енергетичний спектр шуму; 
Sp (и, v) = |F(u, о)? - енергетичний спектр оригінального зображення; 

С(и, v) - Фур'є - перетворення спотвореного зображення. 


Причому останній вираз у (8) появляється у наслідок того, що добуток комплексного 
числа на Його комплексно-спряжене рівне квадрату модуля. Фільтр, що представлений у 
середні дужок часто називають фільтром мінімального середньоквадратичного відхилення або 
вінеровським фільтром. 

Як видно з першої стрічки рівняння (6), проблема нулів у спектрі спотворюючої функції 
при використанні вінерівського фільтру не виникає, за виключенням тих випадків коли 
функції Н(и, v) 1 5, (и, v) перетворюються у деяких точках в нуль одночасно. 

Відновлення зображення у просторовій області отримується шляхом зворотного 
перетворення Фур'є від функції Ё(и, v). Відмітимо, що якщо шум рівний нулю, то і його 
енергетичний спектр перетворюється в нуль i вінерівська фільтрація зводиться до інверсного 
випадку. 


6. Результати досліджень. 


У представленій роботі розглянуто нанесення гаусівських шумів на зображення а також 
відновлення зображень із використанням середньогеометричного фільтру, фільтру середньої 
точки та фільтру Вінера 

Результат фільтрації зашумленого початкового зображення, здійсненого за допомогою 
виразу представлено на рис.5. Розмір використаної маски становив 3x3 пікселі. Як видно з 
цього рисунку 5a, використання середньогеометричного фільтру призводить до ефекту 
згладжування, що порівняльне з дією середньоарифметичного фільтру [1], однак при цьому 
губиться менше число дрібних деталей. Так видно, що верхня частина зображення стала більш 
однорідною і гладкою, хоча при цьому ще можна розрізнити певні елементи поверхні. Крім 
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того, фільтрація приводить до зміни форми гістограми по відношенню як до початкового так 
і спотвореного зображень (Puc3r,46,56). Причому змінюється як форма так і положення 
гістограми. Вона стає вужчою у порівнянні з гістограмою спотвореного зображення в околі 
нульового математичного очікування. Крім того, наній появляється незначний додатковий 
максимум в околі 0,75 


E= 15000 


10000 


5000 F 


а 6 
Рисунок 5. Результат використання середньогеометричного фільтру до спотвореного 
зображення. (а)-відновлене зображення, (6)- його гістограма. 


Виходячи з отриманого результату, за гістограмою можна оцінити ефективність дії 
середньогеометричного фільтру. Для цього скористаємось емпіричним виразом (9): 


0—04—0; , (9) 
тут с, - СКВ спотвореного зображення, а оу - відновленого зображення. 


Аналізуючи дані що передсавлені на рисунках можна показати що під дією 
середньогеометричного фільтру СКВ спотвореного зображення зменшується до величини 0,3 
так що вередньоквадратична величина також зменшується, як показано виразом (10): 


6 10-08 = 0,2. (10) 


Реалізація фільтра серединної точки здійснювалось із використанням виразу (6). Результати 
представлено на рис .6. 


а 6 
Рисунок 6. Результат використання фільтру «серединна точка» до спотвореного зображення. 
(а) -відновлене зображення, (6) - його гістограма. 


Як і у попередньому випадку на відновленому зображенні спостерігається часткове 
покращення. Однак загальна картина залишається зашумленою. Об'єктивну картину змін, які 
наступили після дії фільтру дає гістограмна характеристика, яка показана на Рис. 66. Tyr як i 
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y попередньому випадку основний пік гістограми дещо змістився вправо а його «хвіст» 
простягнувся в сторону великих яскравостей. Тим не менше, порівнюючи СКВ цього 
результату з спотвореним зображенням, можна побачити, що ширина основного піку дещо 
зменшилась , що підтверджується розрахунками виразом (11): 


6 = 1,0 — 0,7 = 0,3. (11) 


Для випадку фільтру Вінера, як робочу формулу для фільтрування використано кінцевий 
вираз (12): 


1 IHGuv)I? 
> Aum бр u,v). а 
H(uv) IH Qu SIC B 


Розрахунок згідно цієї формули є спрощений, оскільки усі величини y ній є наперед 
відомими: H (u, v) - Фур'є образ гаусівського нормального шуму, а Sp (м, v) = |N (u, v)|* - його 
енергетичний спектр; 5/(и, v) = |F(u,v)|* - енергетичний спектр оригінального зображення; 
С(и, v) - Фур'є-перетворення спотвореного зображення. Тому, після зворотного перетворення 
Фур'є, слід очкувати доброго відтворення. Дійсно, як видно з рис. 76 відновлене зображення 
практично співпадає з оригінальним, що підтверджується і наведеною гістограмою. 
Визначення різниці СКВ між зашумленим i відновленим зображеннями становить, як показано 
виразом (13): 


6 = 1,0 – 0,57 = 0,43. (13) 


Це говорить про те що фільтр Вінера у випадку заданих параметрів дає хороші результати. 
Разом з тим, як 1 у попередніх випадках відновлення приводить до появи «хвоста » у правій 
частині гістограми, однак він по інтенсивності набагато менший і створює незначну «засвітку» 
на відновленому зображенні. 


10“ Гистограма зображення після фільтру 


0.5 


а 


Рисунок 7. Результат використання фільтра Вінера до спотвореного зображення. (а) - 
відновлене зображення, (6) - його гістограма. 


7. Перспективи подальшого розвитку досліджень. 


Як видно з отриманих результатів, у методі просторового відновлення, де 
використовується фільтр - маска, результат відновлення зображення є досить незначним 
навіть у випадку спотворення найпростішим шумом, яким є гаусівський. Так для 
середньогеометричного фільтру величина параметру К у відновленому зображенні становить 
лише на 0,2, а для фільтру «серединна точка» - на 0,3. Крім того, у відновленому зображенні в 
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обох випадках має місце «хвіст» на гістограмах в області високоенергетичних значень 
пікселів. Тим не менше візуально спостерігається певне згладжування зображення. Тобто ці 
методи можуть використовуватись у тому випадку , коли спостерігача не цікавлять дрібні 
деталі, а звертається увага на великі області. 

Більш складними є частотні методи, до яких відноситься розглянута у роботі фільтрація 
Вінера. Цей метод добре працює у тому випадку, коли відомі параметри шуму та 
оригінального зображення. Так у нашому випадку величина параметру R = 0,43. При цьому 
спостерігається досить добре відновлення зображення. 

Однак спектр неспотвореного зображення рідко буває відомим. У цих випадках, коли 
спектри шуму 1 неспотвореного зображення невідомі і не можуть бути оцінені, 
використовується підхід, у якому використовується наближений вираз (14) [14,15]: 


^ - 1 |H(u,v)|2 
F(u, v) Е С т) с (и, v) Й 19 


де К - певна константа. Цю константу підбирають емпірично, виходячи з мети 
відновлення зображення. 


8. Висновки 


У даній роботі аналіз проводиться із використанням високорівневої мови програмування 
Python з його бібліотеками для роботи з зображеннями Pillow та OpenCV. запрограмовані 
основні методи відновлення зображень у просторовій та частотній областях. Незважаючи на 
те, що ці методи відрізняються математичними підходами, їх ефективність можна порівняти 
шляхом побудови гістограм та вимірювання їх параметрів. При цьому, можна використати 
параметр R, який у даному випадку, показуватиме відносне звуження гістограми відновленого 
зображення по відношенню до спотвореного, наведено виразом (15): 


в=(1- 2) = 100(o; — оу)/оһ , (15) 


де позначення Ti ж самі що і у формулі (9): о, - СКВ гістограми спотвореного, а оу - 
відновленого зображення. 
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Abstract. When processing and transmitting information in the form of images, the problem of 
reducing their distortion due to various noises is important. Noise reduces the quality of the image 
and thus the perception of the information contained in it. Because of these problems, the ability to 
evaluate the information that can be obtained as a result of analysis using both visual and computer 
methods is impaired. The process of reducing noise in images 15 handled by an area of image 
processing called restoration. Despite the intersection of this area with image enhancement, it should 
be noted that the latter is a largely subjective procedure, while the restoration process is largely 
objective in nature. During restoration, an attempt is made to reconstruct or reproduce the distorted 
image, using a priori information about the phenomenon that caused its deterioration. Restoration 
methods are based on the modeling of distortion processes and the use of reverse procedures to restore 
the original image. In this paper, a new method for evaluating the restoration of distorted images, 
implemented with the help of spatial and frequency filters, is proposed. It consists in converting 
distorted images before filtering into gray scale and obtaining their histograms. The obtained 
histogram is extrapolated by a Gaussian curve and the value of the root mean square deviation is 
determined. A similar procedure is carried out for the restored image. The analysis was performed in 
the Python programming language using the Pillow and OpenCV image libraries. . The restoration 
evaluation parameter R is proposed, which is the ratio of the difference between the root mean square 
deviations of the distorted and restored images to the root mean square deviation of the distorted 
image. The R parameter was estimated for the geometric mean filter, the midpoint filter, and the 
Wiener filter and was found to be equal to 0.2; 0.3 and 0.43, respectively, which correlates with visual 
observations. 
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